Electron temperature in electrically isolated Si double quantum dots 
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Charge-based quantum computation can be attained through reliable control of single electrons in lead-less 
quantum systems. Single-charge transitions in electrically-isolated double quantum dots (DQD) realised in 
phosphorus-doped silicon can be detected via capacitively coupled single-electron tunnelling devices. By 
means of time-resolved measurements of the detector's conductance, we investigate the dots' occupancy 
statistics in temperature. We observe a significant reduction of the effective electron temperature in the DQD 
as compared to the temperature in the detector's leads. This sets promises to make isolated DQDs suitable 
platforms for long-coherence quantum computation. 



Quantum mechanical charge and spin states of elec- 
trons confined in semiconductor double quantum dots 
(DQD) 1-3 have recently attracted much interest, as they 
can be exploited to implement solid-state quantum com- 
putation. One key requirement to perform quantum 
logic operations is a long coherence time for the qubit- 
embodying states. Among other system materials, sil- 
icon is particularly suited to retain spin-coherence for 
long time mainly due to the existence of a stable iso- 
tope ( 28 Si) without nuclear magnetic moment. 4 Another 
approach to mitigate the decoherence introduced by the 
interaction with the environment is the suppression of ex- 
change processes with electrons in the reservoirs. 5 This 
has proven to be beneficial if a chargc-qubit implemen- 
tation is to be preferred. Indeed, single-qubit operations 
have been successfully implemented in trench-isolated sil- 
icon double quantum dots, 6 where electrons are confined 
in a lead-less system and only capacitively coupled to 
gates. 

Here, we investigate the electronic occupation of 
phosphorus-doped isolated DQDs using a single-electron 
tunnelling device (SET) as a charge sensor. By means of 
time-resolved charge detection, the probability of occu- 
pation of each dot is evaluated as the control gate is swept 
across a degeneracy point (i.e. alignment of energy levels 
which leads to electron tunnelling) . From these measure- 
ments we have extracted the electron temperature in the 
isolated quantum system and compared it with the one 
measured in the SET leads. The advantageous effects of 
a lead-decoupled system are quantitatively observed as a 
significant temperature reduction in the DQD. 

The top inset of Fig. 1(a) shows a scanning electron mi- 
crograph (SEM) image of a representative device. The 
use of silicon-on-insulator wafers ensures spatial confine- 
ment of electrons in the vertical direction; reactive ion 
etching is used to define deep trenches which produce 
confinement within the horizontal plane. The silicon 
active layer is doped with a concentration of phospho- 
rus atoms of about 3xl0 19 cm -3 to provide the nano- 
structure with free carriers. From the density of im- 
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planted donors, we estimate that a 40 nm diameter dot 
contains about 2000 electrons. Both electron-beam and 
optical lithography are employed to define the device pat- 
tern. Full details of the fabrication process, tunnel barri- 
ers formation and influence of interface trapping for this 
system are reported elsewhere. 7 In order to keep the op- 
crating point of the detector in the most sensitive region 
while inducing transitions in the DQD, gates G and Gqd 
are simultaneously swept with a fixed ratio. By contrast, 
gate Go is always kept at a fixed voltage. Experiments 
have been performed in a 3 He cryostat at the base tem- 
perature Tfc = 295 mK, unless otherwise stated. 

Fig. 1(a) shows the positions of the detector's Coulomb 
peak maxima as a function of the voltage applied to gates 
Gdd and G. The plot reveals the occurrence of Coulomb 
Blockade (CB) in the SET; the resonances' slope indi- 
cates that both gates are coupled to the SET island 
and significantly affect its potential. In particular, it is 
possible to determine the relative strength with which 
each gate is coupled to the SET in terms of gate ca- 
pacitance ratio. Indeed, the slope of the CB features 
(s = ^Xvcf ' = directly reflects the ratio between 
the SET gate capacitances (see left inset of Fig. 3), which 
can be written as Cq = s- Cq dd \ this defines the voltage 
compensation ratio needed to keep the SET at a fixed 
bias point. Most interestingly, we observe that the tra- 
jectory of the conductance peaks in [Vg,Vg dd ] space is 
nearly periodically shifted. In the bottom inset, the shift 
highlighted by the dashed rectangle is shown as a mis- 
alignment between two adjacent groups of peaks in two- 
gate voltage space. As previously reported, 8 this effect is 
the clear signature of single-charge tunnelling within the 
nearby DQD structure; in particular, for more negative 
Vgdd we expect electrons to move from the top dot (T) to 
the bottom one (B). Each of these events produces a shift 
of about 1% of a CB period which is equivalent to an in- 
duced charge in the SET island of O.Ole, where e is the el- 
ementary charge. This directly depends on the amount of 
mismatch in the coupling of each dot to the detector, as it 
has been found 8 that A oc Cts/C^ t — Cbs/C^ b , where 
A is the modification of the SET island potential upon a 
tunneling event in the DQD, Cts(Cbs) is the coupling 
capacitance between the top (bottom) dot and the SET 
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FIG. 1. (a) Positions of the detector's Coulomb peak max- 
ima as a function of Vg and Vg dd (dotted line). Charge 
transitions in the DQD are sensed as shifts in the position 
of the peak maxima that periodically break their trajecto- 
ries. Dashed lines highlight the loci of the shifts as parallel 
straight lines that separate regions of equal charge distribu- 
tion in the DQD. Top inset: false colour SEM micrograph of 
a representative device with scale bar of 200 nm. Designa- 
tion of electrodes: detector's source (S), drain (D) and gate 
(G); DQD's gates (G DD , G ). Bottom inset: SET current 
in the region highlighted by the dashed rectangle in the main 
plot. Current range: 0pA< Ids <55pA. Arrows highlight the 
peaks' shift, (b) Differential conductance of the detector as a 
function of Vg dd - Vg is simultaneously compensated to keep 
the conductance level at a fixed value. Pairs of arrows high- 
light sharp spikes which reveal sudden change of conductance 
due to tunnelling events in the DQD. 



island and Cs r (Cs B ) is the total self-capacitance of the 
top (bottom) dot. It is worth noting that the shifts have 
all nearly equal amplitude (1%) and draw parallel lines 
(see dashed lines in the main plot) that divide the voltage 
space in regions of stable charge distribution within the 
DQD. This symmetry indicates that the transitions orig- 
inate all from the same spatial location (the DQD) which 
has a fixed capacitive coupling to both the charge sensor 
and the gate electrodes. Fig. 1(b) depicts the differential 
conductance of the detector as the gate voltages Vg dd 
and Vg are mutually compensated along the trajectory 
followed by the lower CB resonance in (a) to keep maxi- 
mum charge sensitivity. Sharp resonances clearly appear 
in correspondence of the voltage shifts; this reflects, once 
more, the occurrence of discrete changes in the electro- 
static potential of the SET due to remote charge transfer. 

In order to observe the time-dependent evolution of the 
DQD near a transition, traces are recorded for different 
gate voltages across a degeneracy point. In Fig. 2 we re- 



port five representative traces of the SET conductance for 
time intervals of 300 s at a sampling frequency of 13 Hz. 
Data are taken in the vicinity of the transition shown in 
the inset of Fig. 1(a). It is evident that near the degener- 
acy point the system switches randomly between the two 
states associated with two discrete conductance levels of 
the SET (i.e. Gt and Gb for an electron occupying the 
top and the bottom dot, respectively). The fraction of 
time that the system spends in each state is a function 
of the gate voltage and approaches 50% at the degener- 
acy point (central trace). Such an information can be 
extracted from the time-dependent plots by counting the 
number of points, n, per conductance bin. A histogram 
of n versus conductance is reported next to each relevant 
trace. The histograms allow one to directly visualize how 
gate voltages control the crossover from localised to de- 
localised charge states. 

The presence of a range of gate voltages for which 
the system is degenerate (as opposed to sharp transi- 
tions at a single-valued voltage) suggests that tunnelling 
is affected by inelastic scattering processes mainly due 
to electron-phonon interactions. The gate voltage spread 
across the transition can be used to extract information 
on the energy level broadening and the effective electron 
temperature in the isolated DQD structure. From his- 
tograms similar to those shown in Fig. 2, the distribution 
functions of Fig. 3(a) can be built by converting the gate 
voltage into energy scale, 
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where chdqd is the lever arm of gate Gdd to the DQD, 
Eq and Vg dd are the degeneracy point electrochemical 



-1 



1 



3 



G R 



«fT I ™ n WW ™ fff 1 f|W Iff W T ™ * 7 llffWfpr^ll ■ 



fLL.il 



G T 
G B 
G T 
G B 

G T 
G B 

G T 
G R 



100 t [s] 200 



300 



FIG. 2. Time-resolved traces of the normalised SET con- 
ductance near a transition point for Vg dd =-2. 549V, -2.551V, 
-2.553V, -2.555V, -2.558V from top to bottom. The conduc- 
tance fluctuates between two discrete levels according to the 
DQD occupancy modification given by the gate voltage ap- 
plied. Histograms on the left hand side account for the time 
the system spends in each state; n is the normalised num- 
ber of points per trace at the conductance level given on the 
vertical axis. 
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potential and voltage, respectively. Following the theory 
by Beenakker, 9 we can treat the discrete energy spec- 
trum in the DQD as a continuum, given that for our 
system Ec ~> ksT^ > Ae, k B being the Boltzmann's 
constant, T®® D the electron temperature in the DQD, 
Ae the single-particle level separation and Eq the dots' 
charging energy. Under this assumption, the distribution 
function can be approximated by the Fermi-Dirac (FD) 
distribution, 
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from which the effective electron temperature can be ex- 
tracted. It is worth pointing out that our system largely 
satisfies the classical limit approximation as Ec ~ 17 
meV, k B T^ D m 30 neV and Ae w 20 ^eV. Fig. 3(a) de- 
picts three representative distributions for different val- 
ues of Th. Solid line fits allow one to estimate the val- 
ues of T®Q D which appear not to be affected by the 
base temperature up to approximately 500 mK, as shown 
in Fig. 3(b) (squared data points). Indeed, it is clear ^ 
from the fitting lines in blue that for T b < 500 mK, E, 
t dqd _ 500 mK w hereas for T b > 500 mK, T b « T° QD . ^ 
This is the evidence that the lowest limit for the electron 
temperature in the DQD is set to roughly 500 mK by 
our cryogenic set-up. In order to correctly extract the 
electron temperature from the occupancy distributions, 
it is paramount to know precisely the value of oidqd to 
be used in Eq. (1). Given that direct transport mea- 
surements are not allowed due to the isolated nature 
of the DQD, the estimate of oldqd has been attained 
as follows. Initially, we have used a value of first at- 
tempt by assuming olbqd — &set — 0.033 eV/V, where 
a set is the lever arm of gate G to the SET. This can 
be easily extracted from the diamond plot of the detec- 
tor. Next, from the observation that T®Q D is nearly 
constant up to T b « 500 mK, we have fitted the occu- 
pancy distribution at T b = 500 mK using cxdqd as a 
fitting parameter and T^"^ D = 500 mK as a fixed pa- 
rameter. The final value of the lever arm so obtained is 
a DQD = 0.055±0.001 eV/V. The fact that the line of best 
fit in Fig. 3(b) for Tf, > 500 mK satisfies the condition 
T b = T®Q D is a decisive consistency check which corrob- 
orates the approach used to evaluate the lever arm. 

It is now interesting to compare the value of mini- 
mum electron temperature for the DQD and the detec- 
tor's leads. To this end, we measure the full width at half 
maximum (FWHM) of the Coulomb peaks in the SET as 
the base temperature of the cryostat is increased. We find 
that the width scales linearly with temperature 10 from a 
minimum value T e det w 650 mK, which is an estimate 
of the electron temperature in the detector (see dotted 
data points in Fig. 3(b)). The discrepancy between the 
two electron temperatures is a convicing manifestation 
of the beneficial effects of decoupling quantum systems 
from leads electrically connected to room temperature 
electronics. However, the presence of other concurrent 
heating mechanisms in the DQD is indicated by the fact 
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FIG. 3. (a) Occupation probability function extracted from 
time-resolved traces as a function of the electrochemical po- 
tential relative to Eq. Data points are taken at different base 
temperatures as reported in the legend. Solid lines are best 
fits according to the FD distribution. Left inset: equivalent 
circuit indicating the main coupling capacitances to the detec- 
tor as well as the self-capacitances of the two dots. Coupling 
to gate Go is omitted, (b) Estimates of electron temperatures 
in the DQD (squares) and SET (dots). Data points for the 
DQD are obtained by using T®® D as a fitting parameter for 
distributions of the type shown in (a). FWHM of the SET 
is evaluated by fitting Coulomb resonances with Lorentzian 
peak functions. Dashed lines indicate the knee points which 
estimate the minimum electron temperatures in the detector, 
r e det « 650 mK, and the DQD, T° QD w 500 mK. 



that the lowest T^ D is larger than any base temper- 
ature value below 500 mK, which is unexpected for an 
isolated system in thermal equilibrium with the phonon 
bath. A number of reasons may be the cause of this ef- 
fect. Firstly, the presence of the charge detector in the 
vicinity of the DQD may produce heating exchange due 
to phonon-mediated back-action. 11 In other words, the 
current flowing through the SET produces ohmic dissi- 
pation in the leads that can be absorbed by the phonon 
bath and transferred to the electrons in the DQD. If this 
mechanism turned out to be dominant, we would ob- 
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serve modifications in the value of T®Q D as a function 
of the detector's bias point. This will be the focus of 
further investigations. Moreover, the weak coupling be- 
tween the electrons and the phonon bath, which is typ- 
ically in place in nano-scale systems because of finite- 
size effects, 12,13 should be considered. This may lead 
to exceedingly long thermalisation times on the scale of 
the time-resolved experiments presented here and conse- 
quently result in an increase of the electron temperature. 
Finally, a temperature difference may also arise due to 
some difficulty in attaining premium thermal conduction 
between the condensed 3 He and the cold finger where the 
device is plugged in. 

In conclusion, this work has shown that control of 
charge state can be achieved in fully electrically isolated 
DQD systems. Most interestingly, we have demonstrated 
that the benefit of decoupling the DQD from the leads 
can be quantitatively estimated as a strong reduction in 
electron temperature. Given that the mitigation of deco- 
herence sources is a crucial requisite for reliable quantum 
manipulations, this system promises to be a viable can- 
didate for charge-based quantum computation. 
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